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Based on recent lattice QCD (LQCD) results obtained at finite temperature, we discuss modeling of the 
hadronic phase of QCD in the framework of hadron resonance gas (HRG) with discrete and continuous mass 
spectra. We focus on fluctuations of conserved charges, and show how a common limiting temperature can 
be used to constrain the Hagedorn exponential mass spectrum in different sectors of quantum number, through 
a matching of HRG and LQCD. For strange baryons, the extracted spectra are found to be consistent with 
all known and expected states listed by the Particle Data Group (PDG). The strange-mesonic sector, however, 
requires additional states in the intermediate mass range beyond that embodied in the database. 

PACS numbers: 12.40.Yx, 12.40.Nn, 14.20.Jn, 14.40.Df 


I. INTRODUCTION 

The thermodynamics of the confined phase of QCD is com¬ 
monly modeled with the hadron resonance gas (HRG) [1- 
8]. The equation of state for strongly interacting matter at 
finite temperature is well described by this model, formu¬ 
lated with a discrete mass spectrum of the experimentally con¬ 
firmed particles and resonances. This finding was verified 
by recent results of lattice QCD (LQCD) [9-12]. However, 
LQCD also reveals that, when considering fluctuations and 
correlations of conserved charges, there are clear limitations 
in the HRG description [11]. This is particularly evident in 
the strange sector, where the second-order correlations with 
the net-baryon number xbs or strangeness fluctuations xss 
are larger in LQCD than those in the HRG model [10, 11]. 
Such deviations were attributed to the missing resonances in 
the Particle Data Group (PDG) database [11]. 

Different extensions of the HRG model have been proposed 
to quantify the LQCD equation of state. They account for 
a possible repulsive interaction among constituents and/or a 
continuously growing exponential mass spectrum [5, 6, 8, 13]. 
The latter was first introduced by Hagedorn [14] within the 
statistical bootstrap model (SBM) [15-17], and was then 
studied in dual string and bag models [18-20]. For large 
masses, the Hagedorn spectrum p{m) is parametrized as 
p(rn) ~ where Th is the Hagedorn limiting tem¬ 

perature and a is a model parameter. 

The main objective of this paper is to analyze LQCD data 
on fluctuations and correlations of conserved charges within 
the HRG model. In particular, we examine whether the miss¬ 
ing resonances contained in the asymptotic Hagedorn mass 
spectrum are sufficient to quantify LQCD results. We focus 
on the susceptibilities xbs and xss^ where LQCD indicates 
the largest deviations from HRG, in spite of their agreement 
on the equation of state in the hadronic phase. 

To calculate fluctuations of conserved charges within HRG, 
one needs to identify the hadron mass spectrum for different 
quantum numbers. For a continuous mass spectrum p{m), this 
issue was addressed in Refs. [21] and [22], where the param¬ 
eters of p{m) in different hadronic sectors were extracted by 


fitting the spectra to the established hadronic states in the PDG 
database [23]. It was shown in Ref. [22] that the Hagedorn 
temperatures for mesons and baryons are different, 
with > r|. The ~ 140 MeV found in [21] is clearly 
below the LQCD crossover temperature = 155(1) (8) MeV 
from hadronic to quark-gluon plasma phase [24—26]. This, 
however, is inconsistent with LQCD, as it implies a large fluc¬ 
tuation of the net-baryon number deep in the hadronic phase, 
which is not observed in lattice simulations. 

In this study we have reanalyzed the Hagedorn mass spec¬ 
trum in different sectors of quantum number, in the context 
of the PDG data, and have shown that there is a common 
Hagedorn temperature for mesons and baryons in different 
strange sectors. We have applied our newly calculated p{m) in 
the HRG model to explore different thermodynamics observ¬ 
ables; in particular, fluctuations of conserved charges. The 
results are compared with LQCD for the strangeness, net- 
baryon number fluctuations, and for baryon-strangeness cor¬ 
relations. We show that HRG, adopting a continuous mass 
spectrum with its parameters fitted to the PDG data, can par¬ 
tially account for the missing resonances needed to quantify 
LQCD results. 

To fully identify the missing resonance states, we motivate 
a matching of LQCD and HRG to extract a continuous mass 
spectrum p{m). In the strange-baryonic sector, this p{m) is 
shown to be consistent with all known and expected states 
listed by the PDG. However, the mass spectrum for strange 
mesons require some additional resonances in the intermedi¬ 
ate mass range beyond those listed in the PDG compilation. 

The paper is organized as follows: In Sec. II, we introduce 
the HRG thermodynamics with a discrete mass spectrum. In 
Sec. Ill, we discuss HRG model comparisons with LQCD. 
In Sec. IV, we extract the continuous p{m) in different sec¬ 
tors of quantum number and discuss fluctuations of conserved 
charges in conjunction with LQCD findings. Finally, Sec. V 
is devoted to summary and conclusions. 
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11. EQUATION OF STATE OF HADRONIC MATTER 


To formulate a phenomenological model of hadronic mat¬ 
ter at finite temperature and density, one needs to identify the 
relevant degrees of freedom and their interactions. In the con¬ 
fined phase of QCD the medium is composed of hadrons and 
resonances. 

The HRG model, in its simplest form, treats the medium 
constituents as point-like and independent [1]. Thus, in such 
a model setup, the interactions of hadrons and the resulting 
widths of resonances are neglected. Hence, the composition 
of the medium and its properties emerge through a discrete 
mass spectrum, 

djS (m - m»), (1) 


where di = (2Ji-|-l)is the spin degeneracy factor of a particle 
i with mass rrii and the sum is taken over all stable particles 
and resonances. 

The mass spectrum in Eq. (1) can be identified experimen¬ 
tally or can be calculated within LQCD. In both cases our 
knowledge is far from complete. LQCD can determine the 
masses of hadronic ground states and low-lying excited states 
with fairly high precision [27]. However, the higher excited 
states are still not well controlled in lattice calculations. 

The spectrum of experimentally established hadrons, sum¬ 
marized by the PDG [23], accounts for all identified par¬ 
ticles and resonances, i.e., confirmed mesons and baryons 
granted with a three- or four-star status, of masses up to 
Tom — 2.4 GeV and ruB — 2.6 GeV respectively. The in¬ 
vestigation of higher excited states remains a significant chal¬ 
lenge for the experiments due to the complicated decay prop¬ 
erties and large widths of the resonances. 

Instead of the hadron mass spectrum (1), the medium com¬ 
position can be characterized by the cumulant [21] 

(m) = diO {m — rrii), ( 2 ) 


such that 


P 


HRG 


dm 


( 3 ) 


Thus, counts the number of degrees of freedom 

with masses below m. 

Since the spectrum (1) is additive in different particle 
species, it can be decomposed into a sum of contributions 
from mesons and baryons, as well as a sum of particles with 
definite strangeness. 

Figure 1 shows the cumulants in different sectors of 
hadronic quantum number with inputs from the PDG. The cu¬ 
mulant of all hadrons is seen in Fig. 1(a) to rapidly increase 
with mass. For m < 2 GeV such increase is almost linear, 
indicating that the hadron mass spectrum is exponential, as 
predicted by Hagedorn in the context of SBM [14, 15]. 

A rapid increase in the number of states is also seen, in 
Figs. 1(b) and 1(c), to appear separately for the mesonic and 


baryonic sector, as well as for the strange and non-strange 
mesons with m < 2 GeV. Baryons of different strangeness, 
as illustrated in Fig. 1(d), follow a similar trend with the ex¬ 
ception of I S'I = 3 baryons, which consists only of fl hyper- 
ons. 

For an uncorrelated gas of particles (and antiparticles) with 
a mass spectrum p{m), the thermodynamic pressure P = 
P/T'^ is obtained as 

PiT,V,p.) = ± J dmp{m) J ^ 

X [ln(l±Ae-'^)-f ln(l±A-ie-"")] , (4) 

where p = p/T, rh = m/T, e = \/fP P mP', and the (±) 
sign refers to fermions and bosons respectively. For a parti¬ 
cle of mass m, carrying baryon number B, strangeness S and 
electric charge Q, the fugacity A reads 

A(r, p) = exp {BpB + Sps + Qpq) , (5) 

where (1 = p/T. Note that, for scalar particles with vacuum 
quantum number, the antiparticle term should be dropped to 
avoid double counting. 

Expanding the logarithm and performing the momentum 
integration in Eq. (4) with the discrete mass spectrum 
in Eq. (1), one obtains 


P = 


E di 

27r2 ^ 
i k—1 


(± 1 ) 


fc+1 


k‘^ 


mf K2{kmi) 


(6) 


where the first sum over i includes the contributions of all 
known hadrons and antihadrons, and K 2 is the modified 
Bessel function. The upper and lower sign is for bosons and 
fermions respectively. The Boltzmann approximation corre¬ 
sponds to retaining only the first term in summation. 

The thermodynamic pressure in Eq. (4), through the mass 
spectrum p, contains all the relevant information about the dis¬ 
tribution of mass and quantum number of the medium. Thus, 
it allows for the study of different thermodynamic observ¬ 
ables, including fluctuations of conserved charges. 

Furthermore, one can turn this argument around, and ex¬ 
plore the implication of the thermodynamic observables on 
medium composition. This approach has been applied, for ex¬ 
ample, in the pure gauge theory to extract an effective glueball 
mass and spectrum based on the lattice results on pressure and 
trace anomaly [28, 29]. For the case of QCD, the constraint 
imposed by the trace anomaly on the hadronic spectrum has 
been investigated [30, 31]. In this work, we focus on the im¬ 
pact of recent LQCD data on the fluctuations of conserved 
charges. 


III. HADRON RESONANCE GAS AND LQCD 

LQCD provides a theoretical framework to calculate the 
equation of state and bulk properties of strongly interacting 
matter at finite temperature. The first comparison of the equa¬ 
tion of state calculated on the lattice with that derived from 
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FIG. 1: (Color online) Cumulants of the PDG mass spectrum in different sectors of quantum number: (a) all hadrons; (b) mesons and baryons; 
(c) mesons of different strangeness; (d) baryons of different strangeness. The lines are obtained from the fit of Eqs. (9) and (12) to the PDG 
data with the parameters listed in Table I (see text). 


Eq. (6) has shown that thermodynamics of hadronic matter is 
well approximated by the HRG with mass spectrum generated 
on the corresponding lattice [2, 3]. 

Presently, we have a rather detailed picture of the ther¬ 
modynamics of hadronic matter, thanks to the advancement 
of LQCD simulations with physical quark masses and to 
the progress in extrapolating observables to the continuum 
limit [12, 32, 33]. Thus, a direct comparison of the equation 
of state from Eq. (4) and LQCD can be performed with the 
physical mass spectrum [4, 10, 34]. 

In Eig. 2(a) we show the temperature dependence of the 
thermodynamic pressure obtained recently in lattice simula¬ 
tions with physical quark masses [12, 32]. The bands of 
the LQCD result indicate the systematic errors due to contin¬ 
uum extrapolation. The vertical band marks the temperature 


Tc = 155(1)(8) MeV, which is the chiral crossover tempera¬ 
ture from the hadronic phase to the quark-gluon plasma [26] . 
These LQCD results are compared in Eig. 2(a) with predic¬ 
tions of the HRG model for the mass spectrum (1), which 
includes all known hadrons and resonances listed by the 
PDG [23]. 

There is a clear coincidence of the HRG and LQCD re¬ 
sults on the equation of state at low temperatures. The pres¬ 
sure is strongly increasing with temperature towards the chiral 
crossover. This behavior is well understood within HRG as 
the consequence of growing contributions from the escalating 
number of higher resonances. 

Although HRG formulated with a discrete mass spectrum 
does not exhibit any critical behavior, it nevertheless repro¬ 
duces remarkably well the lattice results in the hadronic phase. 
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FIG. 2: (Color online) Lattice QCD results of HotQCD [10, 32] and Budapest-Wuppertal Collaborations [9, 12] for different observables in 
dimensionless units: (a) the thermodynamic pressure; (b) the net-baryon number fluctuations Xbb. Also shown are the HRG results for the 
discrete PDG mass spectrum (dashed line) and for the effective mass spectrum in Eq. (11), which contains a continuous part to describe the 
effects of massive resonances (continuous line). 


This agreement has now been extended to the fluctuations and 
correlations of conserved charges [9, 10, 35, 36]. 

In a thermal medium, the second-order fluctuations and cor¬ 
relations of conserved charges are quantified by the general¬ 
ized susceptibilities 

\xy Q» r, ^ 5 ( ' ) 

where {x, y) are conserved charges, which in the following 
are restricted to the baryon number B and strangeness S. 

For HRG with a discrete mass spectrum of Boltzmann par¬ 
ticles, Xxy is obtained from Eq. (6) as 

= dimjK2 {nii) x,y,. ( 8 ) 

Ax=A!/=o ^ 

The susceptibilities (7) and in particular (8) are observables 
sensitive to the quantum numbers of medium constituents. 
Thus, x^y can be used to identify contributions of different 
particle species to QCD thermodynamics [35, 36]. 

Recent LQCD calculations of the HotQCD Collabora¬ 
tion [10] and Budapest-Wuppertal Collaboration [9, 12] pro¬ 
vide results on different fluctuations and correlations of con¬ 
served charges. Thus, the apparent agreement of HRG and 
LQCD, seen on the level of the equation of state, can be fur¬ 
ther tested within different hadronic sectors [10]. 

In Figs. 2(b) and 3 we show the LQCD results on the 
fluctuation of net-baryon number, strangeness, as well as the 
baryon-strangeness correlations. They are compared to the 
HRG model, formulated with the PDG mass spectrum. From 
Fig. 2(b), it is clear that the net-baryon number fluctuation 
in the hadronic phase are well described by HRG, whereas 


strangeness fluctuation xss in Fig- 3(b) and xbs correlation 
in Fig. 3(a) are underestimated in the low-temperature phase. 

Following an analysis of the relations between differ¬ 
ent susceptibilities of conserved charges, it was argued in 
Ref [10] that deviations seen in Fig. 3(a) can be attributed to 
the missing resonances in the strange-baryonic sector. In view 
of Fig. 3(b), similar conclusion can be drawn for the strange 
mesons. 

In general, the contributions of heavy resonances in HRG 
are suppressed due to the Boltzmann factor. However, the rel¬ 
ative importance of these states depends on the observable. 
In the hadronic phase, the pressure is dominated by the low- 
lying particles. At temperature T = 150 MeV, the contribu¬ 
tion to the pressure from particles and resonances with mass 
M > 1.5 GeV is of the order of 7%. However, in the fluctua¬ 
tions of the net-baryon number and baryon-strangeness corre¬ 
lations, such a contribution is already significant and amounts 
to 26% and 33%, respectively. 

Contributions from missing heavy states could be the po¬ 
tential origin of the observed differences between LQCD re¬ 
sults and HRG predictions on fluctuations and correlations of 
conserved charges in the strange sector, shown in Figs. 3(a) 
and 3(b). 


IV. HAGEDORN MASS SPECTRUM AND LQCD 
FLUCTUATIONS 

A. The Hagedorn mass spectrum 

To account for the unknown resonance states at large 
masses we adopt the continuous Hagedorn mass spectrum 


-HRG 

A-xy 
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FIG. 3; (Color online) As in Fig. 2, but for baryon-strangeness correlations xbs (a) and for strangeness fluctuations xss (b). Also shown are 
the corresponding results obtained from the least-square fits to lattice data up to T ~ 156 MeV. 


with the parametrization 

p^{m) = 


Oo 


(m^ 


, 2 \ 5/4 


^m/Tn 


and its corresponding cumulant 



( 9 ) 


( 10 ) 


where is the Hagedorn limiting temperature, whereas oq 
and Too are additional free parameters. 

In general, the parameters of p{m) can be calculated within 
a model, e.g., in SBM [15, 17]. In the following, we constrain 
the Hagedorn temperature and the weight parameters (oq , mo) 
in Eq. (9) based on the mass spectrum of the PDG and the lat¬ 
tice data. In addition, we assume that the same exponential 
functional form holds separately for hadrons in different sec¬ 
tors of quantum number, i.e., for mesons and baryons with or 
without strangeness. 

The analysis of the experimental hadron spectrum, in the 
context of Hagedorn exponential form, has been extensively 
discussed in the literature [13, 15, 16, 37]. In one of the re¬ 
cent studies [21, 22], it was shown that, in htting the Hagedorn 
spectrum to experimental data, one arrives at different limit¬ 
ing temperatures for mesons, baryons, and hadrons with dif¬ 
ferent electric charges. In particular, with p{m) from Eq. (9), 
the limiting temperature for mesons, ~ 195 MeV, was 
extracted to be considerably larger than that for baryons, 
~ 140 MeV. Such Hagedorn limiting temperatures, how¬ 
ever, are inconsistent with recent lattice results, which show 
that the changes from the hadronic to the quarks and gluons 
degrees of freedom in different sectors appear in the same 
narrow temperature range of the chiral crossover. Thus, the 
Hagedorn temperature of baryons should appear beyond the 


chiral crossover, i.e., T§ > 155 MeV. In addition, the LQCD 
data on xbb are consistent with the discrete PDG baryon mass 
spectrum up to T ~ 160 MeV. This seems to suggest that 
large contributions form heavy resonances are not expected in 
Xbb atT| < 155 MeV. 

Prom the above one concludes that it is very unlikely for the 
Hagedorn limiting temperatures in various hadronic sectors to 
differ substantially. Moreover, they are expected to be larger 
than the chiral crossover temperature. Consequently, the ex¬ 
tracted Hagedorn temperature ~ 140 MeV for baryons in 
Refs. [21, 22], though mathematically correct, is disfavored 
by LQCD. 

The reason for the very different Hagedorn temperatures for 
mesons and baryons is that the extraction of the parameters 
in Eq. (9) has been performed over the whole mass range of 
the PDG data. The low-lying baryons drives the fit towards a 
lower Th, resulting in the deviation of Hagedorn temperatures 
among different sectors. 

To avoid the above problem, we adopt Hagedorn’s idea to 
treat the contributions of ground state particles' separately 
from the exponential mass spectrum. In addition, we start the 
continuous part of the spectrum from the onset of the hrst res¬ 
onance in the corresponding sector. Therefore, we apply the 
following mass spectrum 


p(rn) = diS{m — nii) -F p^(m)9{m — m^), (11) 

i 


^ Particles that do not decay under the strong interaction. In this context, 
there are no ground states in the |S'| = 2 and |5| = 3 baryonic sectors. 
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Fit to PDG 

Fit to LQCD 


mo [GeV] 

oo(mo) [GeV^/'^] 

mo [GeV] 

ao(mo) [GeV^/"] 

pH 

0.529(22) 

0.744(40) 

0.425(24) 

0.573(36) 

Pb 

0.145(23) 

0.135(7) 

0.078(13) 

0.108(9) 

pr° 

0.053(8) 

0.064(12) 



pi^~^ 

0.051(12) 

0.046(6) 

0.193(96)(122) 

0.067(27) 

p%^~^ 

1.453(441) 

0.023(20) 

2.469(456)(297) 

0.091(47) 

p%-^ 

0.00194(0) 

0.00027(0) 



Pm 

0.244(17) 

0.341(19) 



pfr° 

0.183(19) 

0.212(17) 



pIT~^ 

0.183(43) 

0.060(9) 

0.378(32)(95) 

0.099(24) 



Constraint 


nix [GeV] 


PH 

0.77526 

18 

pB 

1.2320 

28 

pr° 

1.2320 

40 

p%^~^ 

1.3828 

20 

p%^~^ 

1.31486 

2 

p%-^ 

1.67245 

4 

Pm 

0.77526 

18 

pfr° 

0.77526 

14 

pfr”" 

0.89166 
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TABLE I: (Left) Parameters of the Hagedorn mass spectra in Eqs. (9) and (11), in different sectors, obtained from fits to PDG and LQCD 
data. The Hagedorn temperature has been set to Th = 180 MeV. Sectors of all hadrons, all mesons, and nonstrange mesons include both the 
particles’ and antiparticles’ contributions. In matching the LQCD results, the data for pressure and second-order fluctuations are compared with 
the HRG model through Eqs. (4) and (15). Also shown are the errors of mo arising from the least-square fits, which induce the uncertainties 
in ao(mo) through Eq. (14). For the sectors , and p^j“ , the systematic errors m the approximation schemes are also included 

(see text). (Right) The constraint on the continuous mass spectrum in each sector, given in Eq. (13). 


and the corresponding cumulant 

m 

N{m) = di9(m — rrn) -|- 9(m — m^) / dm (m), (12) 

where {m) is given by Eq. (9). The index i counts the 
hadronic ground states, i.e., states with masses less than rux 
of the first resonance in the corresponding channel. 

With such a prescription for analyzing the hadronic spec¬ 
trum, we find no practical advantage in treating the continuous 
(to) as a two-parameter function of (toq, uq). We therefore 
impose the following constraint on the continuous mass spec¬ 
trum 


N^^^{mx) = Jdmp^im). (13) 

0 

This way, p^ is reduced to a function of a single parameter 
Too . The parameter oq can be determined by 


ao{mo) = A^™‘^(TOa;) 



^m/TH 


(14) 


The above spectrum can now be compared with the experi¬ 
mental data listed by the PDG, in different sectors of quantum 
number. From the analysis of the mass spectrum parameters 
of all hadrons, we find that the best description is obtained 
with Th — 180 MeV. This value is consistent with that re¬ 
cently found in Ref. [13]. In addition, Th — 180 MeV is the 
largest temperature obtained as the solution of the Bootstrap 
equation [38]. In the following, we apply the same Th for 
strange and nonstrange hadrons. 

In Fig. 1 we show that the spectra of PDG hadrons in dif¬ 
ferent sectors are indeed properly described by the asymptotic 


mass spectrum (11) with a common Hagedorn temperature 
Th — 180 MeV. The weight parameters (ao,TOo) in Eq. (9) 
are determined by the composition and decay properties of 
the resonances, hence, they are distinct for each hadronic 
quantum number. The optimal sets of parameters of p(rn) 
in Eq. (9) are summarized in Table 1. The corresponding mass 
spectra are shown in Fig. 1 as continuous lines, whereas cir¬ 
cles indicate the lowest masses rrix of the corresponding fit. 
Also shown, as broken lines in Fig. 1, are the extrapolated 
cumulants below rUx ■ 


It is important for the decomposition of the hadron mass 
spectrum (11) into different sectors, using parameters from 
Table I, to produce results that are thermodynamically consis¬ 
tent. Thus, e.g., the total pressure obtained from Eq. (4) 
with the mass spectrum from Fig. 1(a) should be consistent 
with the sum of meson Pm and baryon Pb pressures, calcu¬ 
lated with the mass spectra in Fig. 1(b). Similar results should 
hold for the pressure when adding up the contributions from 
strange particles in different sectors. This consistency check 
provides further constraints on the mass spectrum parameters 
presented in Table I. 


With the PDG mass spectrum extrapolated to the contin¬ 
uum, we can now test whether heavy resonances can reduce 
or eliminate the discrepancies between HRG and LQCD on 
baryon-strangeness correlations and strangeness fluctuations, 
seen in Figs. 3(a) and 3(b). 


The second order cumulants Xxy, at vanishing chemical po- 
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FIG. 4: (Color online) Cumulants of the discrete PDG mass spectrum (black dashed line) and the corresponding fits (red dashed line) for (a) 
strange baryons and (b) strange mesons. Also shown are the cumulants containing in addition the unconfirmed states (broken-dashed line). 
Continuous lines are obtained by matching the LQCD results to the continuous mass spectra through Eq. (15), assuming that the missing 
strange baryons come solely from the l^l = 1 sector [scheme (I)] or jS'j =2 sector [scheme (II)] (see text). 


tential, are obtained in HRG as 


OO 

Xbb = J ^ pB{m)m^K2 (m) , 


(15a) 


Xss 


OO 

/ 


dm 


Pm + '“(m) 


w? K 2 {rh), 


-H 

Xbs — — 


OO 

/ 


dm 




in K 2 (rh), 


(15b) 

(15c) 


using the mass spectrum p{m) in Eq. (11) and the parameters 
presented in Table I. 

In Figs. 2 and 3, we show the contribution of the contin¬ 
uous Hagedorn mass spectrum to the pressure and different 
charge susceptibilities. The difference between the full line 
(fit to PDG) and the dashed line (PDG) comes from the inclu¬ 
sion of heavy resonances. The results in Figs. 3(a) and 3(b) 
indicate that heavy resonances can capture, to a large extent, 
the differences between HRG and LQCD for strangeness and 
baryon-strangeness correlations. However, at low tempera¬ 
tures, Xbs still differs from the lattice. These deviations sug¬ 
gest that there are additional missing resonances in the PDG 
data in the mass range m < 2 GeV, as they begin to contribute 
substantially to xbs xss ^1 lower temperatures. 


B. Spectra of strange hadrons from LQCD 

To identify the missing strange resonances in the Hagedorn 
mass spectmm, we use the LQCD susceptibility data as input 
for Eq. (15) to constrain p{m) in different sectors. 


We begin with the strange-baryonic sector. The xbs data 
alone do not allow for a unique determination of the contribu¬ 
tion from a particular sector. This is because the observable 
depends only on a linear combination of the spectra, namely 
Pb ~ Pb^~^ + + ^Pb^~^- principle, this problem 

could be resolved with additional lattice data on higher-order 
strangeness fluctuation, e.g., xbbss ™d the kurtosis. 

For our purpose of analyzing the present data, we instead 
make the following simplification. The [S'] = 3 sector is re¬ 
stricted to those states listed by the PDG. We then make the 
assumption that the additional strange baryons come solely 
from the [S’] = 1 [hereafter named scheme (I)] or |5'| = 2 sec¬ 
tor [scheme (II)], and treat the remaining one with the spec¬ 
trum fitted to the PDG.^ The resulting spectrum parameters 
for both schemes are presented in Table I. 

In Fig. 4(a) we show the cumulants of the lattice-induced 
p{m) under both schemes, together with the experimental 
spectra including the unconfirmed states from the PDG. The 
mass spectrum extracted with scheme (I) is seen in Fig. 4(a) 
to follow the trend of the unconfirmed states of the PDG. This 
is not the case for scheme (II). Hence, the extra PDG data for 
the unconfirmed hyperons support the former scenario. 

The Xss fluctuation in general receives contributions from 
both the strange mesons and the strange baryons. However, 
due to Boltzmann suppression, we expect the observable to be 
dominated by the mesonic contribution. This can be inferred 
from the fact that xss ^ Xbs- To be definite, we fix the 
strange baryon contribution to xss by the scenario dictated by 


^ The errors induced by the use of the PDG-fitted spectra are introduced as 
systematic errors for the spectrum parameters. See Table I. 
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scheme (I). The lattice data on xss then directly determines 
the strange-mesonic spectrum. The resulting parameter is pre¬ 
sented in Table 1, with the corresponding cumulant shown in 
Fig. 4(b). 

Similar to the case of strange baryons, the spectrum de¬ 
termined from LQCD requires additional states beyond the 
known strange mesons. From Fig. 4(b), we find that it ex¬ 
ceeds even the trend set by the inclusion of the unconfirmed 
resonances. This may point to the existence of some uncharted 
strange mesons in the intermediate mass range. 

In addition, the general conclusion of an enhanced 
lattice-motivated strange spectra, relative to the PDG, does 
not depend on the chosen functional form of the continuous 
spectrum (9). This follows from the observation that lattice 
data show a stronger interaction strength in the strange sector 
than that expected from a free gas of known hadrons. Within 
the framework of HRG this implies an increase in the corre¬ 
sponding particle content. 

Nevertheless, it is important to bear in mind that such con¬ 
clusion, based on an ideal resonance-formation treatment of 
the hadron gas, is not definitive. For example, the contribution 
to xss from the non-strange sector is also possible through the 
vacuum fluctuation of ss mesons. Such an effect is neglected 
in the current model and a theoretical investigation is under 
way. 

V. SUMMARY AND CONCLUSIONS 

Modeling the hadronic phase of QCD by the hadron reso¬ 
nance gas (HRG), we have examined the contribution of heavy 
resonances, through the exponential Hagedorn mass spectrum 
p{m) ~ to the fluctuation of conserved charges. A 

quantitative comparison between model predictions and lat¬ 
tice QCD (LQCD) calculations is made, with a special fo¬ 
cus on strangeness fluctuations and baryon-strangeness corre¬ 
lations. 

We have reanalyzed the mass spectrum of all known 
hadrons and resonances listed in the Particle Data Group 
(PDG) database. A common Hagedorn temperature, Th — 
180 MeV, is employed to describe hadron mass spectra 
in different sectors of quantum number. This value of 
Th exceeds the LQCD chiral crossover temperature = 
155(1)(8) MeV. The latter signifies the conversion of the 
hadronic medium into a quark-gluon plasma. 

Applying the continuum-extended mass spectrum calcu¬ 


lated from the PDG data, we have shown that the Hagedorn 
asymptotic states can partly remove the disparities with lat¬ 
tice results in the strange sector. 

To fully identify the missing hadronic states, we perform 
a matching of LQCD data on strangeness fluctuations and 
baryon-strangeness correlations with HRG. The parameters of 
the Hagedorn mass spectrum p{m) are well constrained by 
LQCD data in different sectors of strange quantum number, 
using the same limiting temperature Th — 180 MeV. 

The mass spectra for strange baryons inferred from the ex¬ 
isting LQCD data are shown to be consistent with the trend 
of the unconfirmed resonances in the PDG. This is not the 
case for the strange-mesonic sector, where the corresponding 
pirn) exceeds the current data of the PDG, even after the un¬ 
confirmed states are included. This may point to the existence 
of some uncharted strange mesons in the intermediate mass 
range. Clearly, new data and further lattice studies are needed 
to clarify these issues. Moreover, such missing resonances 
could be important for modeling particle production yields in 
heavy ion collisions. 

It would be interesting to assess the effects of resonance 
width on the Hagedorn spectrum. Recent studies suggest that 
the implementation of low-lying broad resonances in thermal 
models must be handled with care [39, 40]. The impact on 
the global spectrum and consequently the thermodynamics is 
currently under investigation. 
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